arXiv:1503.01597vl [cond-mat.mes-hall] 5 Mar 2015 


Lasing in circuit quantum electrodynamics with strong noise 

M. Marthaler^, Y. Utsumi^ , D. S. Golubev^’^ 

^Institut fiir Theoretische Festkorperphysik, Karlsruhe Institute for Technology, D-76128 Karlsruhe, Germany 
^Department of Physics Engineering, Faculty of Engineering, Mie University, Japan 
^Institut fur Nanotechnologie, Karlsruhe Institute of Technology,D-76021 Karlsruhe, Germany 
‘^O.V. Lounasmaa Laboratory, Aalto University School of Science, Einnland 

We study a model which can describe a superconducting single electron transistor (SSET) or a 
double quantum dot coupled to transmission-line oscillator. In both cases the degree of freedom is 
given by a charged particle, which couples strongly to the electromagnetic environment or phonons. 

We consider the case where a lasing condition is established and study the dependence of the average 
photon number in the resonator on the spectral function of the electromagnetic environment. We 
focus on three important cases: a strongly coupled environment with a small cut-off frequency, a 
structured environment peaked at a specific frequency and 1/f-noise. We find that the electromag¬ 
netic environment can have a substantial impact on the photon creation. Resonance peaks are in 
general broadened and additional resonances can appear. 


I. INTRODUCTION 

In recent experimentsii^ an effective single-atom maser 
was realized using a superconducting single electron tran¬ 
sistor (SSET), which is essentially a charge qubit with 
an additional applied transport voltage. The two qubit 
states are coupled to a transmission-line resonator and 
a third state plays a role in the pumping cycle [Fig. 
[1] (a)]. This system has been theoretically studied^^— 
and it has been experimentally shownii^ that lasing can 
be achieved. Theoretically it has also been shown that 
a SSET can be used to create a non-classical photon 
distribution in a resonatoriSiS. The same design can 
be used with a mechanical oscillatorsiMi instead of a 
transmission-line resonator and has been considered as 
one of the possibilities to observe non-classical states in 
such a macroscopic objeciji^. If a Josephson junction is 
directly shunted by a resonator it is also possible to ob¬ 
serve enhanced photon creationi^, multistabilityi^, and 
non-classical effectsi^^— . 

A very similar situation has been studied using a 
double quantum dot coupled to a transmission-line 
resonator-iSd^ [Fig. [T] (b)]. Here an applied transport 
voltage drives a current through a double dot, where 
the energy difference between left and right dot can 
be used to create photons. The correlation of photon 
emission and transport properties has been extensively 
studieci22ri2^. It has also been theoretically shown that 
lasing can be achieved32i2i. Already, in experiments a 
linewidth narrowing has been observed^^ and even clear 
experimental evidence for lasing has been shown^S. Many 
theoretical studies have been performed for this setup. It 
has e.g. been used as a toy model for the quantum pho¬ 
tovoltaic effec1j21 and it has been considered as a hybrid 
microwave-cavity heat engine^^. 

In both cases, SSET and double dot, we consider 
the charge degree of freedom. Charge is strongly cou¬ 
pled to any perturbation from the environment. The 
charge variable of the SSET is strongly coupled to fluc¬ 
tuations in the background charge and the electromag¬ 


netic environmen t^^ 1 ^° . In the double dot, the electron 
is strongly coupled to phonons^. In most theoretical 
studies the effects of noise are modeled by using a lowest 
order master equation, which assumes a smooth spec¬ 
tral function and relatively weak coupling to the envi¬ 
ronment. However, it has already been shown that if 
strong coupling to the bath and higher order transition 
processes are considered, new effects like enhanced pho¬ 
ton emission^^ and inversionless lasing^ can appear. 

We are going to consider strong coupling to the en¬ 
vironment with a structured spectral function. The 
approach used in this paper is similar to the descrip¬ 
tion of strong dephasing in a quantum dot coupled to 
a cavity^. However we consider a system where noise 
creates both dephasing and relaxation, depending on the 
bias point. Additionally we consider charge tunneling 
(either quasiparticles for the SSET of electrons for the 
double dot), which serves as pumping mechanism for the 
three level maser. We will focus on three specific ex¬ 
amples of noise spectral densities: low frequency noise 
created by a strong external resistor—, finite frequency 
resonances2£i21 and 1/f noise^^— . Our results suggest 
that low frequency noise was an important factor in the 
realization of the single artificial atom maser using an 
SSET- and that spurious resonances can lead to addi¬ 
tional resonance peaks. In certain limiting cases it is 
possible to describe the system by a Markovian master 
equation, that treats photon creation as stochastic pro¬ 
cess. Despite the fact that in this case all quantum coher¬ 
ence is lost, the classical coherence of radiation, which is 
the hallmark of lasing is still preserved and we will explic¬ 
itly show the connection of our approach to a standard 
description of the laser— in section [V] 

In the next section we will explain in detail the model 
we use to describe the SSET coupled to the resonator 
and to strong charge noise. We write the Hamiltonian 
in the basis of the charge qubit eigenstates and use a 
standard description of a bosonic reservoir, coupled to 
the island charge. Then we discuss the consequences of 
strong charge noise and introduce the diagonalisation of 
the longitudinal charge noise component. Noise creates 


2 



FIG. 1. A schematic depiction of the lasing cycle for the SSET 
and for the double dot. (a) A cooper pair tunnels trough the 
left junction onto the island and via two consecutive single 
charge tunneling events into the right lead. The energy for 
the process is provided by the transport voltage V, and the 
gate voltage Va can be used to change the gate charge SNq. 
The energy of the island with two additional charges is tuned 
to be smaller than the state with zero additional charges. The 
energy difference, is approximately in resonance with a cavity 
coupled to the island, (b) An electron tunnels from the left 
lead into the left dot. Then the electron tunnels from left to 
right dot, emits a photon and finally tunnels into the right 
lead. 


transport voltage V, the current through the device. The 
Josephson coupling Ej of the junctions should be weak 
compared to the charging energy scale, Eq = I2Cy. 
(Ce is the total capacitance of the island), and the su¬ 
perconducting energy gap A. It leads to coherent Cooper 
pair tunneling, with pronounced consequences when two 
charge states differing by one Cooper pair are nearly de¬ 
generate. In addition, quasiparticles tunnel incoherently 
(with rate oc V/eR^ where R is the normal state resis¬ 
tance on the junction) when the energy difference be¬ 
tween initial and final state is sufficient to create a quasi¬ 
particle excitation, i.e., when it exceeds twice the gap 
(assumed equal for electrodes and island), |Ai?| > 2A^. 
The gate charge Nq = CaVale is chosen to be close to 
one, 5Ng = Nq — 1 < 1/2. This means that we tune the 
SSET close to the JQP cycle where current is transported 
in the following way: A Cooper-Pair tunnels through the 
left junction onto the island, increasing the island charge 
A by 2. This is followed by two quasiparticle tunneling 
events through the right junction, that brings the island 
back to its original configuration^ [Fig. [l](a)]. 

To create a laser-like situation in the SSET we uti¬ 
lize two charge states |A = 0) and \N = 2). These two 
states are coupled by Josephson tunneling through the 
left junction. The similarities to quantum optical sys¬ 
tems becomes more obvious in the base of the charge 
qubit, where we diagonalize the coupled charge states. 
Additionally it is necessary to create population inver¬ 
sion. This is achieved by using a third state, the odd 
charge state \N = 1). 

We describe the system by the following Hamiltonian, 

R — R-sys T Renv E Rcip E Rdiss 

+ Pcoupl + kqp + Pdiss- ( 1 ) 


transitions between the states of our system. We will dis¬ 
cuss these rates and introduce a master equation in the 
Lindblad form which governs the dynamics of our system. 
In a short section we will than show how exactly the same 
model can be used to describe a double dot coupled to a 
resonator. The next two sections will be used to describe 
the effects of coupling to a strong external resistor, to a 
single mode in equilibrium and to 1// noise. In the end 
we summarize our results and we shortly discuss their 
consequences in regard to experiments on mechanical os¬ 
cillators coupled to single charge devices and for further 
work on lasing. 


The system is given by a qubit and a single mode oscilla¬ 
tor in the electromagnetic resonator, Rsys = Aqbit+Aosc- 
The first contribution, Aqbit, describes an isolated qubit 
in the form 

2 

Aqbit = ^ Ec{Rg — N)'^Pn,n 

N=0 

-^{P2,oEPo,2), ( 2 ) 

where we introduced projection operators, 

Pn,m = \N){M\ (3) 


II. MODEL HAMILTONIAN 

In this section we will first discuss the Hamiltonian of 
a SSET in great detail. The SSET consists of two super¬ 
conducting leads coupled by tunnel junctions to a super¬ 
conducting island. A gate voltage Vq shifts the electro¬ 
static energy of the island and controls, together with the 


and \N) are the charge states of the island. 

The two charge state N = 0 and N = 2 are coupled 
by Josephson tunneling through the left junction with 
strength Ej [Fig. [T](a)]. Across the right junction we ap¬ 
ply a large voltage drop, therefore the Josephson tunnel¬ 
ing for this Junction can be neglected^. We diagonalize 
the qubit part of the Hamiltonian Aqbit by introducing 
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the states, 

It) =cos||7V = 0) + sin||7V = 2), (4) 

It) = sin||7V = 0) - cos||7V = 2), 

The states 11) and 11) play the essential role to gener¬ 
ate lasing and \N = 1) is needed to create population 
inversion. It is convenient to introduce Pauli matrices, 


with inverse electron temperature j3e- The density of 
states of quasiparticles is gapped. 


P(w) = 

i 


^i,r) — PO 


0{1 - AVa;2) 

v'l- A2/w2 ’ 


( 12 ) 


and thus only for eV > 2A quasiparticle tunneling and 
simultaneous change in the charge states occurs. In the 
qubit basis, it can be rewritten as. 


V — T 

rqp — -tqp 


COS 2 + -^1,4.) + 


^. = it)(ti-i;)ai, (5) 

= crq- -b (T_, 
cr+ = crl = |t)(-H . 

and the diagonal Hamiltonian of the qubit in this basis; 

i?qbit = ^-^-ScPl,l, (6) 

^E= ^E‘] + 4.{Ec6NgY, 

where we shifted the origin of energy. The Hamiltonian 
of the electromagnetic resonator is 

Hose = Wo Ya . (7) 


+ sin-(P4.,i -b Pi,^) 


H.c. 


(13) 


For 5Ng > 0, the condition cos(^/2) > sin(^/2) is sat¬ 
isfied and the probability for the pumping process | I 
,n) —>■ |l,n) —>■ |t,u) [Fig. [2Kb) left panel] is higher than 
that for the anti-pumping process Itj'ra) —>■ |l,u) —>■ |4,,n) 
[Fig. [2Kb) right panel]. At the symmetry point, SNq = 0, 
both process occur equally cos(^/2) = sin(^/2) = l/-\/2 
and thus no population inversion is created. 

The qubit is capacitively coupled with the 
transmission-line oscillator. It is also coupled with 
environmental charge fluctuations. These couplings 
appear in the same form as. 


The eigenstate of the system Hamiltonian Hsys is char¬ 
acterized by the qubit states, a =t 5 4.j !> and the number 
of photons n as |cr, n) with the energy E^ + nuJo where 
P^/ 4 , = ±AEj2 and Ei = —Ec- For lasing, the qubit is 
operated around the resonant condition, 

duj = 0, i5a; = Wo — AE, (8) 

so that the states | 4-, u -b 1) and | f, n) are degenerate 
[Fig. [2] (a)] and thus coherent emission or absorption of 
photons can be induced. 

For photon emission, population inversion has to be 
created by quasiparticle tunneling^i^. The Hamiltonian 
of the quasiparticles reads 

~ (9) 

r—RJ i 


Koupi = [g (a'^' -b a) -b x] (Po.o - ^ 2 , 2 ), (14) 

i 

where we introduce another oscillator bath to model the 
environmental charge fluctuations, 

77env = (15) 

i 

The bosonic bath is characterized by the correlation func¬ 
tion, 

(x(t)x(O)) = ly ^A(w)n-(w)e—‘. (16) 

The equilibrium distribution function for bosons reads. 


where the subscript R denotes the right lead, while the 
subscript I denotes the superconducting island. We only 
apply a large voltage drop across the right junction, 
therefore we neglect quasiparticle tunneling through the 
left junction^. The effect of the tunneling of quasiparti¬ 
cles on the qubit is described by the following Hamilto¬ 
nian 


Iqp — Tqp {P2,l + T’l.o) + H.c., 
i,k 

The quasiparticles obey the Fermi distribution 


( 10 ) 


n (w) 


1 

1 - e-^"^ ’ 


(17) 


where /3 = l/ksT and T is the noise temperature. The 
spectral density. 


A(w) = 27r^|tjp[d(w-Wi)-(5(w-bWi)], (18) 

i 


can be arbitrary. This model is versatile and able to 
describe a variety of charge noise sources in equilibrium. 
For example, if the source of the environmental charge 
fluctuation is the external circuit the spectral function is 
given as. 


( 11 ) 


A(w) = wReZ{(jj)/R k ■ 


(19) 
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where Z(oj) is the effective impedance as seen from the 
quantum system^i^, and Bx is the resistance quantum. 
The electromagnetic environment is known to have a sig¬ 
nificant impact for single charge devices^. This model 
is also able to describe the coupling to another mode in 
the electromagnetic resonator with frequency ujl, which 
is off-resonant and thus in equilibrium. The spectral den¬ 
sity is, 

fV(w) = a;Lec[(5(a; - wl) - i5(a;-I- wl)] , (20) 

where ec characterizes the coupling strength. For an ex¬ 
tensive discussion of the form of the effective impedance 
and its relation to the coupling to the electromagnetic 
environment see Ref. [45l| . 

In the basis of the diagonal qubit Hamiltonian, the 
couplings between the qubit and the oscillator consists 
of a longitudinal and a transverse part, Hcoupi = Vz + V±. 
The longitudinal coupling reads, 

Vz = [g [a) + a) + x] cos^a^, (21) 

and we separate the transverse coupling into two parts 
V± = Hch + Vg with 

Hch = a: sin^CTa;, (22) 

Vg ~ g(a^(7- + acr+) sin^, 

where we performed the rotating wave approximation, 
namely we neglect terms aV+ and aa- connecting the 
states I t) n) and | |, n — 1) energetically separated by 
AE + LJo- Now Vg induces transitions between the qubit 
states I t,'ra) and | 4,,n -|- 1) [Fig. [5] (a)]. The transverse 
coupling to the bosonic bath 14h causes relaxation [Fig. [5] 
(a)], while the longitudinal coupling 14 cause dephasing. 
The relaxation effect 14h is weak since the states | f, n) 
and I i, n) are energetically separated by the gap AE and 
the relevant frequencies of the spectral densities consid¬ 
ered in this paper are all smaller than AE. Further by 
tuning to the symmetry point SNq = 0 (^ = tt/2), the 
longitudinal coupling Vz disappears. However, for lasing, 
in order to create population inversion, the qubit has to 
be tuned away from symmetry, SNq > 0. 

Finally, we also account for the mixing of the single 
mode in the resonator to an external circuit again mod¬ 
eled by oscillators, 

Hdiss = J2^f4d^, (23) 

i 

Vdiss = + H.c.. 

i 

Since the precise form of the spectrum is not relevant 
we simply assume it smooth around the frequency of the 
resonator, which leads to a decay rate of the form 

«: = y'dt^|tfp(d.(t)4(0))e-‘>‘, (24) 

i 

The temperature in the oscillator reservoir does not have 
to coincide with the noise temperature and since it does 
not substantially effect our results we choose the temper¬ 
ature in the oscillator to be small. 



5Ng 5Ng 


FIG. 2. (a) Energy diagram for the SSET and the elec¬ 

tromagnetic resonator with n (left panel) and n -|- 1 photons 
(right panel) at the resonance condition 5uj = u)o — AE = 0. 
Two states | 4) ’i') and | 4, n -f 1) are energetically degener¬ 
ate. The odd charge states |l,n) and |l,n -|- 1) are shown 
at the energy —Ec + eV. The longitudinal coupling to the 
bosonic bath Fz, which causes dephasing, increases away from 
the symmetry point SNq = 0. It may spoil the trans¬ 
verse coupling to the oscillator Vg crucial for coherent las¬ 
ing. The impact of the longitudinal coupling to the bosonic 
bath Vch, which causes relaxation, will be small since the 
states I 4i S'lid | s'-re energetically separated by the 

gap AE. (b) Creation of population inversion by quasiparti¬ 
cle tunneling. For SNq > 0, the condition cos(5/2) > sin(^/2) 
is satisfied. The probability for the pumping up process 
I 4, fi) I ti fi) (ffift panel) is higher than that for 

the pumping down process | ti|li^) I 45 ^) (right 
panel). At the symmetry point, 5No = 0, both processes 
occur equally, cos(^/2) = sin(^/2) = and thus no pop¬ 

ulation inversion is created. 


A. Lasing-like behavior in incoherent photon 
emission 


In order to generate lasing, first the qubit is tuned to 
the resonant condition, where j f, n) and j 4, -I- I) are 

energetically degenerate [Fig. [5] (a)], and then photon 
emission or absorption is possible. This condition must 
be fulfilled at a bias point 5Nq > 0, such that population 
inversion is created. At the same time, it inevitably cre¬ 
ates longitudinal coupling between the qubit and charge 
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noise T4 ■ The size of the fluctuations is characterized by 

cos^^ {x{tf) = J ^ N{uj)n-{uj). (25) 

Then the longitudinal coupling with the bosonic field 14 
competes with the transverse coupling with the oscillator 
Vg. We consider the case when the fluctuations in the 
level splitting can be larger than the Rabi frequency, the 
transverse coupling between oscillator and SSET, 


5\/Wsin^ < \/(a^(0^)cos^, (26) 


where (n) = (a^a) is the average number of photons in 
the system. As we tune away from degeneracy the trans¬ 
verse coupling to the oscillator that is needed to create 
photons, decays as sin^ at the same time, the fluctuations 
in the level splitting grows as cos 4 We diagonalize the 
longitudinal part 14, using a polaron transformatio n"^^4^ , 

17 = exp[-jcos^(p-l-p')(T2]. (27) 


Here, in order to write the resulting Hamiltonian in a 
compact form we introduce the operators. 




4 {bj - h) 


2 u)i 


/ . 5 (a 

p = i 


t _ 


a) 


UJQ 


(28) 


After the transformation the coupling between SSET the 
oscillator and charge noise takes the form 

Ecoupi = C/14oupiC/t = 14h + (29) 

where I 4 and 14h are now dressed by polarons as, 

Vg^g sin a V+ -f H.c., (30) 

i4h = sin 5u V+ -f H.c., 

where S = pcos^ and S' = p'cos^. Here 14 may be 
treated as the perturbation. 14h may be also treated 
as perturbation because the relevant frequencies of the 
spectral functions considered are smaller than the energy 
splitting AE. The quasiparticle tunneling is dressed with 
polarons as well and the transformed Hamiltonian reads, 


Eqp = Tqp 


cos|(-e-4«+S')/2p^^^ 


+Pi,te 


i(S+S')/2 


H.c. 


(31) 


However, quasi-particle tunneling processes are domi¬ 
nated by the energy scale of the transport voltage eV 
and phase fluctuations cause only minor effects as we 
will discuss later. 

Here at this point we summarize the non-perturbative 
part Hq and the perturbative part H', 

Hq = Hsys + Henv + T^qp + T^diss > (32) 

H' = Vg-\- 14h + T ^iss ■ 


Then the reduced density matrix is obtained by tracing 
over all reservoir degrees of freedom, consisting of the 
charge noise, quasiparticles and dissipation in the oscil¬ 
lator T/^env T -Pqp T Tf^iss j 


p{t) = TrR 
po oc e~^^° 




(33) 


where the initial density matrix is normalized as Trpo = 

1 . 

Under the situation we are considering, because of 
strong environmental noise, the coherent coupling be¬ 
tween the states 11, n) and 11, n+1) is lost and transitions 
among different photon states becomes stochastic. The 
quantities of interest are the average number of photons 
(n) = (a^a) and the Fano-Eactor 


F = 


(n^) — (n)^ 

(n) 


(34) 


The Fano-Eactor defines the width of the photon number 
distribution. In a pumped oscillator one would expect a 
photon number distribution, similar to a Poisson distri¬ 
bution, in this case the Fano-Eactor should be close to 
F = 1. If we only heat the oscillator and we have a ther¬ 
mal photon number distribution the Fano-Eactor should 
be close to F = (n) -I- 1. 


B. Transition rates 

Let us take a closer look at few selected transition 
rates. We do this to make the connection between our mi¬ 
croscopic Hamiltonian and our Markovian Master equa¬ 
tion clear. First, let us consider the quasiparticle tun¬ 
neling processes for generating the population inversion. 
The transition rate of the process |l,n) —> | 4-1''^) within 
Fermi’s golden rule in I4p reads, 

£ /■“ 

rq,„>^|i,n)=2cos2 -Rey^ 7qp"(i) C'”,i/2 

xCg(t)44edSr,"-Ei,n-eV)t^ ( 35 ) 

where the particle-hole propagator is given by, 

7qp”(w)=(Tqp(t)rqy = J J duj'duJp{uj') 

xpiuj + a;')/-(a; + uj')f+ 

(36) 

where /~ = 1—and /'*' is the Fermi distribution. The 
matrix element is given by C”;^ = |(n|e“*''’'® \n)\^ r; 1 up 
to the order of (p/wo)^- For all further purposes we will 
ignore the effect of higher orders of (p/wq)^ on the quasi¬ 
particle tunneling rate, since it has no significant impact 
in our regime. For a discussion of effects this coupling to 
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the oscillator can have see e.g. Ref. The correlation 
function Cg(t) is well known from P {E)-Theoiy^; 




(37) 


= exp 


cos^ X 


, N{uj) 

duj — 5 — 


X ([1 — cos(a;t)] coth(/3a;/2) — i sin(a;t))] , 


It causes a broadening in the energy dependence of the 
tunneling rate by ^ {x{t)^). However, we will always 

operate in the regime eH — 2 A:^cos^ ^/{x(t)^, therefore 
the broadening effect can be neglected. Further, since the 
bias voltage is much larger than the energy difference, 
which is of the order of the charging energy Eq, eV ^ 
~ Ei,n\, the transition rate is approximately, 

r| 4 .,n)^|l,n) ~ COS^(^/2) Fqp , rqp=Gqpy, (38) 

where Gqp = 27r|t'^PppQ is the quasiparticle tunnel con¬ 
ductance. 

The crucial part in our theory is the photon absorp¬ 
tion and emission rates induced by Vg. Although these 
rates are governed by the charge fluctuations a:, it turns 
out that an incoherent pumping by quasiparticle tunnel¬ 
ing also creates additional dephasing and modifies these 
rates considerably. Therefore, we consider second order 
perturbation theory in gsin^, but infinite order in pow¬ 
ers of the quasiparticle tunneling rate Tqp. The photon 
absorption rate is given by, 

nOO 

rq,„).^q,„+i)=25^ sin^^ Re/ dt Cg{t) 

Jo 

x*^l'I;ii)(<)- (39) 

The matrix element reads /"■'■^ = |(n-|-l|e*‘^ a^e®'® |n)p ~ 
n -I- 1 up to the order of (g/ojo)^. From Eq. (l39l) . we can 
see that the correlator Cg{t) causes a broadening of the 
resonance of SSET and oscillator. Therefore we will get 
direct information about the photon emission from the 
correlator. 

The propagator describes the time evolution 

of the off-diagonal component of the reduced density 
matrix (t,n|p| i,n + 1), i.e. the dephasing. Eor ze¬ 
roth order, it simply oscillates with frequency Suj as 
(t) = —One should keep in mind that in 
our system Fqp is the strength of an incoherent pumping, 
that creates an additional dephasing. We have to con¬ 
sider the broadening effects by this additional dephasing 
when the correlator Gg(t) decays weakly or does not de¬ 
cay as is the case when we later relate the environment 
with another mode in an electromagnetic resonator. The 
dephasing can be treated conveniently using the real-time 
diagrammatic technique^^— . The propagator, diagram- 
matically expressed as fig. [3] (a), is written in the Eourier 
space as. 


|4„n-|-l> 


(w) = 


1 


uj + Suj — a, 


|4,,ti-|-1) 


(w) 


(40) 


t 0 


|t.n) 


(a) 


n 


U,n+1 


+ a 


+ a a + 


(b) 


t 0 


a 


bn) |t,n) 



|t,n+l) |bn+l) 


FIG. 3. (a) The propagator describing the time evolution 

of the off-diagonal component of the reduced density matrix 
(t, u|p| 4, u +1). (b) The leading diagrams for the self-energy. 
The two diagrams are second order expansions in Vijp. 


The self-energy can be expanded again in powers of Vg, 
Vch and also I^p. We account for the leading terms of 
the expansion in VJjp, corresponding to two diagrams in 
fig. |3] (b), which actually cause the dephasing^; 

xe*('^‘^+A^/2)*6i(<) -h ^(gVwo) ■ (41) 

As we discussed when we derive Eq. (1351) . the charge noise 
does not influence quasiparticle tunneling for large bias 
voltage eV. Then the self-energy can be approximated 
as ImCTj|’"^j^^(—(5w) ~ —rqp/ 2 , which is independent of 
the gate charge SNq- By neglecting the real part, which 
gives the renormalization of the frequency Suj, we obtain 
the approximate form, 

(42) 


within this approximation, we obtain, 

rIt,n>.^q,«+ 1 >~(?^ + 1) Sg{Suj) = 2 (n + 1 ) g'^ sin^^ 

pOO 

X / d^Gg(^)e*'^‘^*-^■^p‘/^ (43) 

Jo 

where the dephasing effects of the incoherent pumping 
can be simply included via an exponential decay with 
the decay rate rqp/ 2 . 

Within the same approximation, the rate of relaxation 
from I t, n) state to | n) state caused by Vh is calcu¬ 
lated as, 


poo 

rq,n)^|t.n>=2 sin^C Re J dt Gch(t) (t), 


(44) 
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where the correlation function reads (Appendix. A), 
=C'g(t)sin^C [(a;(t)a;(0))-cos2^(p(t)a;(0))^] , 

(45) 


the time evolution in orders of the coupling to the reser¬ 
voirs. Dissipation and Quasiparticle tunneling can be 
considered to be in the Born-Markov limit, therefore we 
truncate the expansion at lowest non-zero order. For the 
charge noise we can truncate at lowest order and perform 
the Markov approximation for 


and 


5g(<5w) < Fqp , S,h{AE) < AE. (49) 


(46) 

Up to the order {g/ujo)"^ the relaxation rate becomes, 

r|4„n)-(-|t.n>~*S'ch(AA) = 2 

pOO 

xRe (47) 

Jo 

The contribution of this term is often small because of 
rapid oscillations with frequency AE, which is of the or¬ 
der of ujQ and larger than any other energy scale but eU. 
Therefore the integral is convergent for t ^ 1/AE. 

The effects of the environment are contained in two 
noise correlators, Cch and Cg [eq. (1571) and (|55]) ]. Both 
correlators depend on the rotation of the qubit. At the 
degeneracy point, ^ = 7r/2, they are reduced to, 

Cg = 1, (48) 

Cch = (x{t)x{0)). 

The first result means that the charge qubit is protected 
against dephasing. The second result means that the 
relaxation from | f) state to | ),) state is caused by the 
charge fluctuations. We will discuss these correlators and 
the corresponding spectral functions, Sg^Suj) and S'ch de¬ 
fined in eq. (H51) and (H7)) respectively, for explicit choices 
of noise spectra later. 

The spectral functions tell us at which energy ab¬ 
sorption (emission) is most effective, however one should 
note that the effect each of these spectra has on the las¬ 
ing properties of the system is markedly different. The 
spectral function Sg{Sijj) describes the dephasing effects, 
which lead to an overall broadening of the resonance be¬ 
tween oscillator and SSET. The corresponding operator 
CT_a^ flips the qubit state, and creates a photon at the 
same time. Therefore the effectiveness of the creation of 
photons is proportional to the size of Sg^Suj). Only for 
large Sg{duj) we will have a significant number of pho¬ 
tons in the oscillator. In contrast to this Sch{AE) has a 
detrimental effect on the creation of photons. The relax¬ 
ation operator tT_ allows transitions between the states 
I t) and II). If therefore S'ch (A A) is large we have an¬ 
other efficient channel, that closes the lasing cycle shown 
in fig. [T] and only few photons will be created. 


C. Master equation 

To find the density matrix in the stationary limit, we 
describe the system by a master equation. We expand 


Then the Master equation reads. 


P — ^ [TIsys; p\ T (Ag T Ach T Aqp Ahiss) P 1 (^1^) 


where the Lindblad operator is separated into 4 parts. 
The coupling between oscillator and SSET is described 
by 


Lgp = 


Sg{6uj) 


( 2 aa+pa-a^ — a^aJaa+p — pa-oJ aa+) 


Sg{—5uj) /„ t + 

H—-(2cr_o'pacr+ — aa+cT-a' p 

—paa+a-oJ) . 


(51) 


Qubit relaxation is contained in the super-operator Lch, 
which reads 


^ch(AE) 

Ach P= -- [2a+pa- - a-a+p - pa-a +) 


Sch(-AE) 


{2a-pa+ — a+G-p — pa+O-). (52) 


Incoherent pumping processes by the quasiparticle tun¬ 
neling is contained in. 


AqpP — (2P,l[_|_j^ 

n=0,l 

~Pn+l,nPn+l^nP ~ PP'n+'i-,nPh+l^n) i (^ 3 ) 
and dissipation in the oscillator part reads, 

AdissP = {2apa) - aJ ap - paJ a) (54) 


Here the Hamiltonian i7sys and the Lindblad operator 
do not mix the off-diagonal and diagonal components of 
the reduced density matrix. Therefore, the system is de¬ 
scribed by the diagonal components and evolves stochas¬ 
tically. In the basis |ct, n) we can than write the equations 
of motion as 


Pi,n — Aqp sin (<|/2)p'f'^^ + Aqp cos (^/2)p^ yj , 

-|-K(n + l)/9i,„+i - (/on + rqp)pi,„ (55) 

P'\,n — Aqp COS (^/2)pi yy -)- Sg(^duj^(^Tl -)- 

-|-S'(-A£’)p4,,„ + K{n + l)/9t-, n+l (Aqp Sin'(e/2) 

+Sg{Su;){n + 1) + S{AE) + kti) , 

Pi.,n — Aqp sin (^/2)pi yj Ag(Jca)np'|'^yy_i 

-|-5'ch(A£')pt.,,y -I- K{n + l)p4„n+i - (r qp cos' (e/2) 
+Sg{du}){n + 1) + S{-AE) + ku) , 

where we used the notation p^r^n = {n,a\p\a.n). As we 
will show in section El we can rewrite this master equa¬ 
tion using standard methods and obtain a form which 







is formally equivalent to the equation of motion for a 
laser—. 

In the next three sections we will discuss the stationary 
solutions of the master equation, first for low frequency 
noise (see eq. (ESD), than for coupling to a spurious reso¬ 
nance (see eq. (|68ll l and at last for 1/f-noise. One should 
note here that we assume that 1 /f-noise can be described 
as Gaussian and ergodic noise, which is not necessarily 
true and depends on the microscopic model^i^. We dis¬ 
cuss the distribution probability of photons, 

the average number of photons (n) = We 

will also investigate the current through the SSET, which 
is given by 

/ = ehqp -b sin^ -b cos^ , (57) 

where the density matrix for the qubit states is obtained 
by tracing out photon number state as, 

Pa =^{(J,n\p\(j,n). (58) 


III. THE DOUBLE-DOT 

In a double dot we consider a situation where a single 
charge can either sit on the left or right dot or both dots 
are unoccupied. Therefore we can write the double-dot 
Hamiltonian as, 

Hqut = elII, 0)(1, 0| + e«|0,1)(0,1| + eo|0, 0)(0,0| 

+t(|l,0)(0,l| + |0,l)(l,0|) (59) 

Here the eigenstates of the double-dot are given by |1,0) 
(left dot occupied), |0,1) (right dot occupied) and |0,0) 
(both dots unoccupied). To map this case to the SSET 
we replace the states of the SSET with the double dot 
states using the following rules, |0,1) —>■ | = 2), |1, 0) —>■ 

|A^ = 0) and |0,0) —>■ |A^ = 1). All energies have to 
be mapped accordingly. The coupling between the dots 
is equivalent to the Josephson coupling, t = Ej/2. A 
detailed discussion of the lasing cwle in the case of the 
couple dot can be found in Refs. ^ andHJ. 

For the double dot it is also reasonable to diagonalize 
the subspace which contains the charge states |1,0) and 
|0,1). This gives us states which act in an equivalent 
fashion to the states defined in eq. o, 

It) =cos||l,0)-bsin||0,1), (60) 

It) =sin||l,0)-cos||0,1), 

tan 4 =-. 

eR — £L 



FIG. 4. The correlator Sch(AE) as a function of the energy 
splitting AE for a fixed noise width y/ksTec = 0.07 (all 
energies are given in units of the charging energy Ec). (blue) 
ec = 0.05, (magenta) ec = 0.03, (brown) ec = 0.01. For the 
parameters ; Ej = 0.2, Fqp = 0.0325. 


This allows us to write the Hamiltonian of the double dot 
in the form, 

AE(j^ 

i7qbit = ^-^+eo|0,0)(0,0|, (61) 

AE= y/4t^ +[cl- en)^, 

The terms in the master equation, which effect the pho¬ 
ton creation (EU and decay dSH) can then be adapted 
directly for the states | f) and |),) of the double dot. The 
pumping rates, which create the population inversion, 
are described by eq. (1531). These terms can be adapted 
for the double dot, using the substitution rules discussed 
above. 


IV. THE ELECTROMAGNETIC 
ENVIRONMENT 

In this section we will consider two examples that have 
been considered previously within the context of tunnel¬ 
ing in small junctions^. The first case is coupling to a 
large ohmic impedance. We have chosen this example 
because it gives results similar to 1/f noise. In fact, as 
we will show in section El one can see 1/f noise as the 
classical limit of a large ohmic impedance. The second 
example is coupling to a single mode in equilibrium. Ef¬ 
fects of this have been measured for a single charge device 
with a well characterized environment e.g. in Ref. Is^ 


A. High Impedance Environment 

If the qubit is coupled to an ohmic resistor the effective 
impedance in eq. (IT^ takes the form^ 

Re^(w) ecujR 

-Rir “ 
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where wn = l/RCes is the RC-cutofT frequency, ec = 
jCeS and CeS is a capacitance that depends on the de¬ 
tails of the coupling to the external resistor with ohmic 
resistance R. We consider noise sharply peaked for low 
frequencies, such that ujr ksT. In this case the de¬ 
phasing correlator becomes, 

cos^ — e e^fc cos^ ^k^Tt 

C^{t) = e .(63) 

The decay time can not be larger than t = l/Tqp. There¬ 
fore we can expand the exponent in this correlator in the 
short time limit for luh <C Tqp and get 

Cg{t) = exp {—ec cos^ ^ [it + ksTt^] ) (64) 

While we will not explicitly consider the opposite limit 
ujR ^ Tqp, one should note that in this case our the¬ 
ory reproduces the standard quantum optics results ex¬ 
pected for a smooth spectral noise density. However, we 
will discuss this transition for small charge noise strength 
ec ^ ^qp’ 

The limit considered here corresponds to a large resis¬ 
tor, i? —>■ oo, in the effective impedance^, 

ReZ{uj)/RK ~ ec6{uj). (65) 




The dephasing spectral function at the energy Soj be¬ 
comes 


Sg(Suj) = 


1 


2 y ec cos^ 


(rqp-iC 5 u; + e^ COS^ 


xerfc 


Tqp - i((5a;-I-ec cos^g) 
2\/ec cos^^fcsT 


-I- c.c. 


( 66 ) 


where erfc[a;] is the complementary error function. We 
can immediately analyze two important limits of this 
spectral function. 


Sg{SLjj) oc 


e 


( 5 u; + e(;; cos^ 5)2 


4 e^ cos 2 5 
2rqp 


cos ^yJecksT » Tqp 
cos ^^ecksT <C Tqp 

(67) 


If we have strong noise we have a Gaussian spectral func¬ 
tion with a width defined by cos^y/ecksT. The Gaus¬ 
sian function is peaked at AE — uj ^ cos^ec- This cor¬ 
responds to the fact, that in each qubit flip, a photon is 
created and additionally the energy ec cos 5 is transfered 
to the environment. Previoulsy it has been shown that in 
a similar situation, such a shift of the resonance can be 
used to create inversionless lasing^ or enhanced photon 
emission in a double doli^. 

Note that the width and center of Sg{duj) depends on 
the rotation of the qubit. Directly at the degeneracy 
point, SNg = 0 , we have cos5 = 0 and therefore Sg is 
simply a Lorentzian with the width being determined by 
Tqp. Overall we see that for strong low frequency noise 


FIG. 5. The average number of photons (n) and the Fano- 
Factor as a function of the energy splitting AE (all energies 
are given in units of the charging energy Ec)- (blue) ec ~ 
0.005, (magenta) ec ~ 0.01, (brown) ec = 0.015. For the 
parameters: Ej = 0.2, ksT — 0.5, uo = 0.4, Fqp = 0.0325, 
K = 0.0006, g = 0.004. 


we get, as one would expect, a Gaussian resonance peak 
and for weak noise we get Lorentzian resonance peak. 

The spectral function Sch{AE) is proportional to the 
relaxation rate. However, for any given configuration of 
our system this correlator will be evaluated for energies 
larger than the Josephson coupling, AE > Ej. Gener¬ 
ally for low frequency noise this means that the correlator 
has already decayed for relevant frequencies and can be 
ignored. However, one should note that we have a sig¬ 
nificant broadening, because of the dephasing caused by 
Tqp, and as we move away from degeneracy, cos 5 > 0, we 
get an additional broadening through cos^^/ecksT and 
an overall frequency shift ec cos 5 (see eq. (IM)l L In fig. |4] 
we can see a plot of the correlator Sc\x{AE). For this plot 
we have chosen a fixed value for the noise width y/ecksT 
and than changed the coupling strength ec- One can see 
that for large coupling the frequency shift has a signif¬ 
icant impact and the spectral function has a peak for 
AE > Ej- However, if the width is mostly caused by the 
noise temperature the correlator decays monotonically. 

That means that the effect of transverse coupling for 
low frequency noise is only important for a very strong 
coupling to the environment. Therefore in our discus¬ 
sion of I/f noise (see sec. |V|, where we make the formal 
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transition to the classical high temperature limit, we can 
ignore transverse coupling because we keep tc and there¬ 
fore the frequency shift small. 

In fig. [5] we show the average number of photons for 
high ohmic impedance, calculated by numerically solv¬ 
ing eq. (1551) in the stationary limit. As predicted from 
eq. (HZl) we see that the resonance peak changes from 
a Lorentzian to a more Gaussian form for increasing 
noise coupling ec- Nonetheless, we still have a signifi¬ 
cant Lorentzian component given by Fqp. One should 
also note that the maximum shifts to energies smaller 
than the oscillator frequency, AA < w. The reason for 
this is the decrease of the coupling strength as we tune 
away from the degeneracy point. We will discuss this 
in more detail in sec. 0 We also show the Fano-Factor 
which shows the standard behavior for a laser. When 
we start to populate the oscillator, the distribution has a 
thermal shape, and therefore the Fano-Factor is large. As 
we move closer to resonance the number of photons grows 
and the distribution takes a form similar to a Poisson dis¬ 
tribution, which means the Fano-Factor moves closer to 
one. 


B. Coupling to a single mode 



AE 



In the limit of coupling to a single mode, the effective 
impedance is reduced to^, 

KeZ{uj)/RK = ec['J(w — ojl) + 5{uj + ujl)] ■ ( 68 ) 

In this case the dephasing spectral function becomes^ 

S'g((5w) = exp (-[ 77 + -b T]-]) (69) 

^^77-77- 2rqp 


n\m\ F^ -(- (Jo; -(- (jT, -(- rn)uJLY ' 


77± = 


ec cos^ ^ 


qp 

±1 


OJL 




This gives us several peaks in the photon number. How¬ 
ever one should note that the largest peak is still lo¬ 
cated at A A — wq = 0. Additionally we have peaks at 
AE—uJo = rhcoL, with fh = ±1, ±2 .... For small temper¬ 
atures, /3 —>■ 00 , the environment can only absorb energy 
and we only have peaks at ttt. > 0 . 

The relaxation spectral function can be be written in 
a compact way by defining 


*S'ch.eff(w) = tan^ ^ 


TTT]- 


rL + (w - ojlY 


(70) 


-b 


TT7J+ 


+ 




^qp + (w -b Wl)^ F^ -b (w -b 2 ujl )‘^ 


V- 


F2 + (a; - 2u:lY 


V-V+ 

Tqp + W" 


If we convolve this spectral function with eq. dUl) , we 


FIG. 6. The correlators Sg{Soj) and Sch{AE) as a function 
of the energy splitting AE, for coupling to a single mode, 
(all energies are given in units of the charging energy Ec)- 
(brown) ec ~ 0.05, (magenta) ec = 0.075, (blue) ec = 0.1. 
For the parameters : Ej — 0.2, ksT = 0.02, Fqp = 0.0325, 
(jjQ = 0.4, ojl — 0.075, g — 0.004. 


obtain the relaxation spectral function 

S'ch(AF;) = exp(-[? 7 +-b 77 -]) (71) 

j.n m 

X 'Y' ~ Sch,es{AE -b (77 -b m)uJL) ■ 

^' n\m\ 

n,m 

Here we can see that the total spectral density S'chi will 
consist of multiple peaks which are repeating with dis¬ 
tance (jJl- This means relaxation will be strong always 
for AE = rhujL, with fh = ±1, ±2,.... Each maximum 
in the relaxation rate will decrease the number of pho¬ 
tons at that particular qubit splitting. At the same time 
incoherent Cooper-pair tunneling is possible and we get 
an additional maximum in the current. Together with 
the maxima which are caused by the lasing transition it 
is possible that we see several peaks in the current. 

In fig. | 6 ] one can see an example for Sg{Suj) and 
Sch{AE) for the case of coupling to a single mode. If 
the system is strongly coupled (ec is large) we get sev¬ 
eral resonance peaks in the coupling between the SSET 
and oscillator. However, these peaks decay very quickly 
as ec becomes smaller. To see clearly visible peaks in 
the relaxation spectral function SchiAE) we would have 
to choose very large ec- This is because we choose a 
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AE 



dephasing correlator in the classical limit 


Cg = exp 



duj S{uj) 


cos(wt) — 

J 


(72) 


where we defined S{uj) = N{u}) coth{u}/2k bT). This is 
equivalent to replacing the coupling to a bosonic reser¬ 
voir in eq. (ITTl) with a coupling to a classical fluctuat¬ 
ing variable. Formally, we can make the transition to 
the classical case by choosing the spectral function to be 
constant, iV(w) = Af, and the temperature to be large. 
In this case the spectral function becomes. 


S{u}) 


ksTAf 

UJ 


(73) 


The factor AfksT characterizes, together with the low 
frequency cutoff Wmin, the strength of 1/f noise. The 
resulting dephasing correlator reads. 


FIG. 7. The average number of photons (n) and the current 
7 as a function of the energy splitting (all energies in units 
of the charging energy Ec). (blue) ec = 0.05, (magenta) 
ec = 0.075, (brown) ec = 0.1. For the parameters: ujl = 0.1, 
Ej = 0.2, g = 0.004, Fqp = 0.0325, k = 0.00005, fcsT = 0.02, 
LJo = 0.4. 


off-resonant mode lol < Ej- 

Similar to the case for low frequency noise, relaxation 
decays quickly for coupling to a single mode. But it can 
still have a significant impact. In fig. [7] we see the av¬ 
erage number of photons (n) and the current / through 
the SSET as a function of the energy splitting. As the 
coupling to the noise ec increases we see the maximum 
number of photons decrease. However, the second reso¬ 
nance peak for AE — 6uj = ujl becomes comparable to 
the first resonance. For large ec the second peak is even 
higher than the first peak. This is the case because the 
contribution of Sch{AE) has decayed. The current shows 
a peak for every peak in the photon number. One can 
also see additional structure in the current which are a 
result of the incoherent cooper-pair tunneling caused by 
Sch{AE). 


V. 1/F NOISE 


The existence of noise with a 1// type spectrum seams 
to be universal to all physical systems and has been an 
object of intense stud y^^'^^i^^~ — . Generally I// noise is 
considered classical noise, therefore we will consider the 


Cg = exp 


/ fceTA/'Reln) 

^minO 4.2^ 

V 2^^ ) 


(74) 


The decay rate of the correlator is much larger then the 
cutoff frequency, this means that the logarithmic time 
dependence ln(a;ininf^) has no impact on the decay of Cg. 
We define a new constant ec = to de¬ 

scribe the noise strength. Here we replaced the time t in 
the In-function with the timescale given by I/Fqp. From 
this we get the dephasing spectral function 


Sg{5u:) = - 


ec cos^ 


(rqp-i(j>u;)2 
_g 4e^ cos^ qT 


(75) 


xerfc 


Fqp - i5uj 
2^Jec cos^ ^knT 


+ c.c. . 


which is the classical limit of eq. (ESI)- 

For classical low frequency noise relaxation will be 
rather small, and decay quickly, as AE increases (see 
sec. IIV Al) . Additionally our dephasing spectral function 
is symmetric, Sg{5uj) = Sg{—Suj). In this case we can 
derive an approximate analytical solution for the density 
matrix in the stationary limit. 

We want to derive an effective equation for the proba¬ 
bility distribution of the number of photons in the oscilla¬ 
tor pn = Pcr.ra- Wc do this by tracing out the degrees 
of freedom of the SSET in the equation of motion given 
by eq. ([31]). Using the relation -|- -|- 

and the assumption that the time scales of the SSET are 
faster than the time scales of the oscillator we can form 
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a closed set of equations , 


dt \ Pi,n ) 


^l,n 


0^2,n 

Pi 


P 2 


f ai,n Sg(6uj)n \ / Pf,n-i \ 
y Sg{Suj)n 02,n J \ Pi,n ) 


I ( PlPn-1 
V P'^Pn 
rqp(cos 2^ + 7) 
4(cos^ + 3) 
rqp(cos2^ + 7) _ 
4(cos^ — 3) 
2rqpCos^(C/2) 

3 + cos ^ 

2rqpSin"(g/2) 

3 — cos ^ 


- Sg{6u})n, 
Sg{Suj)n , 


(76) 


This set of equations can be solved in the stationary limit 
and we get an equation for the effect of the artificial atom 
on the oscillator 


Pn = r+p„_i - (r+_pi+ r„ + Ku) p„ (77) 

+ (r„+l + K{n + 1)) Pn+l , 

where r+ = Ft,™ cos'^(^/2) is the rate increasing the 
number of photons, T" = sin‘*(^/2) is the rate de¬ 

creasing the number of photons and 

Tt.k = 8rqpS'g(^a;)n/ [rqp(cos^ -I- 7) -I- 2ASg{Suj)n]. 

(78) 

For cos(^/2) > sin(^/2) we have a net increase of the 
number of photons and the rate is directly propor¬ 
tional to Sg{6u}). One should note here that eq. (iTTll is 
formally equivalent to the equation of motion as it can 
be derived for the same system in the case of coherent 
coupling using standard methods^. In fact if we choose 
the limit of small charge noise, ec —j' 0, eq. ([77| is ac- 
tually exactly equivalent to the standard lasing result of 

Ref. 0 

Eq. dZll) can be solved and we get for the density 
matrix 

Pn = PoIIj^Q . , (79) 

Jd +1 

_ Fqp cos^(^/2) 

3k 

_ Fqp sin (^/2) Fqp (cos 2^ -|- 7) 

3k 24Sg{Su}) ’ 

where po is a normalization constant. 

If cos(^/2) > sin(^/2) and Sg is not to small this func¬ 
tion has a sharp peak. In this case the average number of 
photons is given in good approximation by the position 
of this peak, 

{n)^ A-B (80) 

_ i^qp e _ Fqp (cos 2^ -I- 7) 

3k ^ 24Sg{5uj) ■ 


We see here that the number of photons increases with 
the size of Sg{5oj), since the second term in eq. (1501) 
decreases. There are two key observations we can make 
from this analytical equation for (n). 

First one should note that in the limit of strong 1/f- 
noise, cos^y/ecksT > Fqp, the resonance peak will be 
shifted to smaller frequencies. This is the case because 
of the decay of the transverse coupling between SSET and 
oscillator with the increase of AE. To find the approxi¬ 
mate position of the actual peak we solve d^ESgiSui) = 0 
for strong 1// noise and get the peak position, 

AEmax ~ ^ fwo + \J^l- l^ecksTj . ( 81 ) 


This shift can also be observed for coupling to a high 
impedance environment [Fig. [5]. 

Now let us consider the number of photons for the 
qubit close to the symmetry point, ^ = tt/2, where we 
assume that it is in resonance with the oscillator, AE = 
Lo. In this case we can expand the number of photons 
around the symmetry point and get, 

(n)^^(^-7r/2) (82) 

Otx 

Eqp / Fqp ksTec 

V 2 Fqp 

We can see here that the second, negative term in this 
equation decreases quadratically, as we move closer to 
symmetry. The positive first term decreases only linearly. 
To create coherent coupling between oscillator and SSET 
it is very important to maximize the number of photons. 
From eq. (I82[) we can immediately see that the number 
of photons has a maximum for 


(C-V2)2) • 


(e-V2) 


g^Fqp 

nksTec 


(83) 


We see here that especially for large charge noise ec it 
can be of advantage to stay closer to the symmetry point. 
The loss in population inversion is compensated by the 
suppression of charge noise. 


VI. CONCLUSION 

In this paper we discuss the single artificial atom maser 
which can be realized using a SSET coupled to an oscil¬ 
lator or a double dot coupled to a transmission-line res¬ 
onator. For such a system it is necessary to bias the sys¬ 
tem away from the symmetry point, to create population 
inversion, which makes it susceptible to charge noise. In 
this case it is a reasonable approach to describe the sys¬ 
tem using a polaron transformation. We have shown the 
transformation of the qubit Hamiltonian and the result¬ 
ing new noise operators. The relevant noise correlators 
have been shown in explicit form and we discuss them for 
two relevant examples, low frequency noise and coupling 
to a single mode. Then we perform the transition to the 
classical case, which describes 1/f-noise. 
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We find that low frequency noise creates a Gaussian 
resonance peak, which is shifted to smaller qubit energy 
splitting. Coupling to an additional mode can create ad¬ 
ditional resonances in the fundamental mode and in the 
current. One should note here that coupling to an addi¬ 
tional mode can even create peaks in the current if there 
is no excitation in the oscillator. 

For our discussion of 1/f-noise we can show that our 
master equation is formally equivalent to the standard 
description of a laser—. We can find an analytical solu¬ 
tion for the density matrix and the average number of 
photons. We can then show that it can be of advantage 
to be closer to the symmetry point, especially if charge 
noise is strong. This has direct implications for experi¬ 
ments on mechanical oscillators coupled to an SSET. To 
create coherent lasing, one should tune the system rather 
close to the symmetry point, since the reduction in pop¬ 
ulation inversion affects the photon number only linearly, 
but the reduction in charge noise is quadratic. 

The model presented in this paper can also describe 
many other system that are currently studied. An ex¬ 
ample are suspended carbon nanotubes. Electrons can 
hop onto a free state of the nanotube, and it has already 
been shown that coupling to the mechanical vibration of 
the nanotube can be used to create cooling^iS. Simi¬ 
larly lasing could be achieved. Since charge is relevant 
degree of freedom for these systems, strong coupling to 
noise should be expected. 
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APPENDIX 


To find an explicit form for the relaxation correlator 
(HSl) we use the generating function. 


where T is the time sorting operator. Using function we 
can show that the relaxation correlator is given by 


Gch(<) 


dPF 


X(t') = (5(t — t' + rj) + 5(t — t' — 77))/2 

+rj) - 5(t' + 77))/2 

= 0 


The explicit form of the generating function reads. 


(85) 


F = exp 



dt2F(ti,t2)&iti - t2) 


( 86 ) 


F = {sit? ^ ^i{ti)^{t2){x{tl)x{t2)) 

+ cos^sin^ fi{ti)X{t 2 ){x{ti)p{t 2 )) 
+ cos^sin^ A(ti)^(t2)(p(fi)i(<2)) 
+ COS^ ^X{tl)X{t2){p{ti)p{t2))) . 


We can apply eq. (1851) to this form of the generating 
function and get 


Cch{t) = sin^ ^ [(a:(t)a:(0)) - cos^ ^{p{t)x{0))'^] (87) 
x(e*PWe-*P(0)). 
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